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A Bose-Einstein condensate (BEC) undergoing measurement-based feedback cooling is analysed 
with a full quantum-field simulation. Two experimental setups are considered: a BEC in a cavity and 
a trapped BEC undergoing phase-contrast imaging. An experimentally important parameter regime 
is found where the simulation results diverge from the predictions of both single-mode quantum 
models and multi-mode semiclassical models. In this regime, simple cooling schemes cannot control 
the quantum noise limited heating due to measurement backaction. We describe a feedback scheme 
that is specifically adapted to the measurement, which can still produce cooling. 
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A Bose-Einstein condensate (BEC) in a dilute atomic 
gas is the premier platform for the investigation of quan- 
tum fields. Since it can be well isolated from its envi- 
ronment and is highly controllable using a combination 
of optical, rf and magnetic fields pQ, it has the poten- 
tial to address a broad range of research questions in 
fundamental and applied science. These include studies 
in quantum nonequilibrium thermodynamics [2], entan- 
glement of massive particles [3] and the quantum simu- 
lation of both cosmological phenomena, such as Hawk- 
ing radiation emitted from a black hole's event horizon 
[3], and phase transitions in condensed matter, includ- 
ing superconductivity and quantum magnetism [5. . Fur- 
thermore, Bose-condensed sources are likely to be key 
components in a range of future technologies, such as 
improved precision inertial sensors based on atom in- 
terferometry, where the sensitivities of current devices 
are limited by the properties of thermal sources [3 [7]. 
This wealth of research opportunity has therefore led 
to great practical interest in the ability to control the 
spatial state of the BEC's quantum field. Research has 
predominantly focussed on open-loop control of the con- 
densate by direct control of optical and magnetic po- 
tentials. In contrast, continuous measurement feedback 
control of BECs is rarely employed due to the perceived 
fragility of the quantum state to measurement backac- 
tion. However, this intuition does not account for the 
advantages active feedback control can provide. Active 
feedback can control properties of the BEC using the 
information gained from the continuous measurement, 
which includes properties that would have been adversely 
affected by the measurement backaction. Furthermore, 
active feedback provides robust and reliable behaviour 
that can not be matched by open-loop control schemes 
[8]. The advantages of feedback controlling a BEC have 
been demonstrated for a single-mode model [TU] and 



for a multi-mode BEC under the semiclassical approxi- 
mation [TTl |T3] . However, an investigation of a feedback- 
controlled Bose gas that incorporates both the full quan- 
tum field and the multi-mode dynamics has not been 
considered. In this letter we show that a correctly de- 
signed feedback-control scheme can indeed cool a degen- 
erate Bose gas to a steady state. Importantly, we show 
that this feedback scheme must be designed with a full 
quantum-field simulation in order to account for the ef- 
fects of spontaneous emission. 

Continuous measurement of BEC has been achieved 
using optical fields, but heating due to spontaneous scat- 
tering of the photons placed strong limits on the life- 
time of the condensate [121 [H] . The fundamental limit 
to this heating can only be improved by placing a BEC 
in a cavity |15H17) . but it is unknown whether using the 
measurement information to perform active feedback can 
overcome this heating and cool a condensate. This prob- 
lem was investigated using Gaussian and semiclassical 
approximations in [TT] [T2 HB] • Under the Hartree-Fock 
semiclassical method, which assumes the condensate has 
a fixed number of atoms, all the same single-particle 
state, it appeared that it was possible to cool a BEC un- 
der the effect of a measurement. However, without a full- 
field calculation it is difficult to determine if the heating 
process has been properly characterised. Full quantum- 
field calculations of BECs undergoing continuous mea- 
surement were made possible by a recent tool based on 
the Number-Phase Wigner (NPW) phase-space distribu- 
tion [TTIH2T] , which allows the simulation of a cooled BEC 
with realistic measurement strengths. 

We examine both a BEC in a cavity and a BEC un- 
dergoing density measurement via phase-contrast imag- 
ing. In both cases, for sufficiently weak measurement 
the Hartree-Fock simulation matches the full-field simu- 
lation and the BEC can be cooled using feedback. How- 
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ever, under more typical experimental parameters, the 
Hartree-Fock simulation fails to predict a heating effect 
caused by the scattering process in the measurement. For 
both systems new controls are proposed, and the full-field 
NPW simulation demonstrates that net cooling can still 
be achieved during high resolution measurement. 

Cavity-mediated Measurement: We consider cooling 
a BEC with a measurement mediated by a cavity, as 
shown in Fig. [T] This feedback mechanism has been 
studied extensively for neutral atoms and optomechan- 
ical resonators [2"2Tf2~i] . Furthermore this measurement 
has been demonstrated experimentally for optomechani- 
cal resonators [25] and ultracold atoms [261 EI] • Light is 
passed through an optical cavity containing the trapped 
BEC. When the light is sufficiently detuned from the 
atoms such that the excited level is minimally populated, 
and the cavity is sufficiently lossy that the cavity dynam- 
ics can be adiabatically eliminated, then the phase of the 
transmitted light contains a signal proportional to the 
overlap of the atomic density with the spatial envelope 
of the optical field in the cavity. A homodyne measure- 
ment of this phase can then be used to continually update 
the best estimate of the atomic state, known as the filter. 
The mean position of the trap can then be moved by an 
amount proportional to the negative of the BEC's bulk 
momentum, which will cool the BEC. The equation of 
motion for the filter that governs the single atom version 
of this system is found in [22] |28] . Generalising this filter 
evolution to apply to a multi-atom quantum field with a 
negligible inter-atomic interactions and converting it to 
Stratonovich form gives: 



d t p = -i{H F (u s ),p] + jV[C s ]p + jC[C 6 ]p+ y/lH[Ct]p V, 

(1) 

where energy is in units of Hu>, position has units of 
Jhj (moj) and time has units of l/u>; p is the condi- 
tional state of the atomic field; 7 is the measurement 
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strength; T>[c]p = cpc 
herence superoperator; Ti\c\p 
the innovations superoperator; C[c]p 



cp + p& 



is the deco- 
- (c + c')/5 is 
c + ct)H[c]p - 

T,H[c 2 ]p+ (&&)p — cp& is the Stratonovich correction su- 
peroperator, and 77 is a Stratonovich differential that is 
derived from the measurement signal in the laboratory. 

= J dx ^ (x)c^(x)^(x) is the measurement operator 
where c^(x) — cos 2 (£2; — f ) defines the intensity of the 
intracavity optical field. The atomic-field Hamiltonian 
Hp = J dx ip* (x)hp(x, u s )4>(x) depends on hp(x,u s ) — 
I + x' 2 ) + u s x(P)/(N), the single-particle Hamil- 
tonian containing the BEC's kinetic energy, potential 
energy due the to harmonic trap, and feedback. P 
I dxt/i' (x)(— iip'(x)) is the full-field momentum operator 
where ip'{x) — d x ip(x), and N — j dx^ (x)ip(x) is the 
multi-particle number operator. £ = 2ttx /X is the ratio 
of the natural length scale of the trap xq and the wave- 
length A of the optical cavity [55] 1 and oj is the angular 
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FIG. 1: Schematic diagram of a BEC under control. In (a) the 
measurement is mediated by a cavity as the homodyne signal 
is proportional to a density moment of the atomic cloud. In 
(b) a direct image of the BEC is taken using phase-contrast 
imaging. In both the controls are applied to the magnetic 
field with the aim of taking the BEC as close as possible to 
the ground state. 



frequency of the harmonic trap. 

In the limit £ <C 1, the cosine-squared term is approxi- 
mately a position measurement. A BEC under a position 
measurement was used as a testing ground for the NPW 
simulation method due to the existence of reliable bench- 
marks [21 . However in contemporary experiments a very 
small £ can be difficult to achieve so we will consider the 
more moderate values of £ = 0.1 and £ = 0.5. We will 
complete simulations using both the Hartree-Fock and 
NPW method. While the NPW method encapsulates a 
more complete description of the physics, the Hartree- 
Fock simulation is faster to perform as it assumes that 
the atoms occupy the same single-particle state. This 
state defines a macroscopic 'wavefunction' an (x) for the 
atoms, which obeys the following evolution equation: 

d t a H (x) = ( - ihp(x) - -/{^(x) 2 - 2c i (x)C i ) 

+V7Q r))a H {x), (2) 

where C\ = J dx c^(x)\aH(x)\ 2 . 

The NPW method provides a full-field solution to the 
master equation given in Eq. (JlJ [21] . The solution con- 
tains a function of the same dimension as the Hartree- 
Fock order parameter, aj^, but physical predictions are 
only obtained by averaging over multiple realisations of 
an extra noise increment 77/. Integrated alongside the 
fields ajv are a set of complex- valued variables w, and 
the lowest order expectation values of operators are ex- 
tracted by the following weighted averages of these func- 
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FIG. 2: (a) Comparison of the Hartree-Fock method versus NPW for simulation of a cooled BEC in a cavity for different £ 
strengths. The NPW representation was integrated using Eq. Q, averaged over 10 paths, and plotted in red with circles. The 
Hartree-Fock representation was integrated using Eq. averaged over 100 paths, and plotted in purple with stars. A small 
value of £ = 0.1 is plotted in (a,i) while a higher strength of £ = 0.5 is plotted in (a,ii). Both plots have 7 = 5, iV = 100 and 
u s = 1. The Hartree-Fock method only agrees with NPW when £ = 0.1, it fails to predict the additional heating present in 
the system when £ = 0.5. (b) A comparison of a controlled BEC under direct measurement without and with the quantum 
noise control averaged over 10 paths. (b,i) shows a comparison of the per particle energy of the system, while the density 
fluctuations in the system without and with control are shown in (b,ii) and (b,iii) respectively. In (b,i): integration of the NPW 
representation Eq. Q without quantum noise control, u c = 0, is plotted with red with circles; integration with the addition of 
the quantum noise control, u c = 5, is plotted in orange with with diamonds. (b,ii) and (b,iii) are contour plots of the BEC's 
density, all lines indicate equidensity parts of the BEC, the colour range of the lines range from (light) blue hues to (dark) red 
hues to indicate low density to high density respectively. All integrations were performed with u s = 1, £ = 0.5, 7 = 5 and 
N = 100. 



tions [20]: 

(a| f) ) - A[aj* } ] (: ha] : sym ) = A[a<a$] (3) 

where :: sym denotes symmetric ordering and A[g(a)] = 
J2n w n9( a n)/ ^ n w n is the weighted average over the 
different stochastic realisations of the noises iff. The 
weighted stochastic differential equations (WSDEs) for 
a BEC in a cavity are [21] : 

d t ottf(x) = -i (h F (x,u s ) + y/jc 6 (x) r)f) a N (x), 

d t w = (-2 7 (Cf-2QA[Q]) + 2V7Q *7)^(4) 

where — J dx c^(x)\aj^(x)\ 2 . 

The comparison of the Hartree-Fock and NPW simu- 
lations [JT] of a BEC in a cavity are shown in Fig. [21 (a), 



where the initial state is a coherent state with a spa- 
tial amplitude given by a displaced Gaussian <p(x) = 
(27r)" 1 / 4 v /iV/crexp {-(x - x ) 2 /4ct 2 ) where x is the 
gaussians displacement from the origin, a is it's vari- 
ance and N is the average total number of the conden- 
sate. When the cosine measurement is close to a position 
measurement, £ = 0.1, we see the NPW and Hartree- 
Fock simulations match well, with both predicting cool- 
ing of the BEC. For larger £ the cosine measurement 
is no longer 'position-like', in that it has some curved 
structure on the scale of the BEC. When £ = 0.5 we see 
the full-field solution predicts significantly more heating 
than the Hartree-Fock solution. If an experiment relied 
on the Hartree-Fock solution they would see unexpected 
irremovable extra heating in their system, demonstrating 
the importance of having these full quantum-field solu- 
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tions when examining quantum control of a BEC. 

The Hartree-Fock method misses the additional heat- 
ing in the BEC because it neglects certain spontaneous 
emission events. The measurement performed on the 
BEC is derived by assuming the atoms that constitute 
the BEC can be described as a two-level system interact- 
ing with off-resonant light. The light is sufficiently off- 
resonant that we can adiabatically eliminate the excited 
state. All the spontaneous emission events in the system 
are now encoded as scattering events between the light 
and BEC. These spontaneous emission events change the 
phase of the outgoing probe beam, which in turn is how 
the cosine-squared moment of the BEC is measured. As 
the Hartree-Fock approximation assumes all atoms are 
in the same single-particle state, it restricts which spon- 
taneous emission events can be described. In contrast, 
the NPW method also allows non-collective states of the 
atoms, and can thereby simulate all spontaneous emis- 
sion events. We will now demonstrate that this new sim- 
ulation method can be used to design and test a control 
scheme that can cool the BEC even in the presence of 
this spontaneous emission noise. 

The fact that the Hartree-Fock simulation worked 
when £ = O.f shows that the net effect of the heating due 
to the spontaneous emission is partially mode-dependent. 
The net effect of the feedback results from the competi- 
tion between the heating due to the scattering and the 
cooling due to the feedback. The feedback in the simula- 
tion was cooling the centre-of-mass motion of the BEC, 
which for £ = 0.1 is very close to the spatial mode of the 
heating due to measurement backaction. When £ = 0.5, 
energy was added to modes other than the position mo- 
ment mode, so we hypothesise that the solution is to add 
an extra control which targets the cosine- squared mode 
of the BEC. 

A methodology used to damp specific modes of a BEC 
is described in [55]. We can apply this method to create 
a control that targets the cosine-squared, c^(x), mode 
of the BEC. The new single-particle Hamiltonian that 
includes this new quantum noise control is 

h FC (x,u s ,u c ) = h F (x,u s ) + c ^ 

x [ dyc^y)lm[(^(y)^(y))}, (5) 



where H F n = Hp + ^ J dx ^ (x) 2 i()(x) 2 contains the 
atomic nonlincaritics, which can dramatically modify the 
behaviour of a trapped BEC. M v (x) = {\i v * ^>^)(x) is 



where u c is the strength of the quantum noise control for 
the cavity-mediated measurement and c'Jx) — d x c^(x). 



NPW simulations of BEC cooling with this new con- 
trol were performed without and with the novel quan- 
tum noise control, and the results are shown in Fig. [2]- 
(b). In part (b,i) we can see the heating due to quantum 
noise is completely cancelled by the quantum-noise con- 
trol. Comparing parts (b,ii) and (b,iii) we see that the 
quantum noise control vastly improves the stability of the 
BEC's density. In (b,ii) we see the quantum noise causes 
the BEC to break apart. In contrast, (b,iii) shows that 
the quantum noise control stops this spread. 



The conditional master Eq. ([T]) that describes the 
cavity-mediated measurement is only valid when inter- 
atomic interactions are negligible. Although a non- 
interacting condensate can be created with a dilute 
atomic sample or via a Feshbach resonance [I . it is 
technically demanding. Furthermore, there are instances 
where large inter-atomic interactions are desirable, such 
as for stable atom laser operation [3jJ] 131] an d for the 
generation of nonclassical atom laser states [35]. For- 
tunately, feedback-control is still possible in this regime 
under a direct measurement of the BEC via off-resonant 
light. 



Direct Measurement : Real-time continuous measure- 
ment of a BEC has already been achieved in the labo- 
ratory using phase-contrast imaging [131 133] . where an 
off- resonant laser beam passes through the BEC, and 
the resulting phase shift is proportional to the atomic 
column density. The conditional master equation for a 
quasi-one-dimensional 'cigar-shaped' BEC under phase- 
contrast imaging is derived in [TT], and in Stratonovich 
form is: 



(6) 

I 

the measurement operator, which is the density oper- 
ator convolved with a kernel function p, v (x) where we 
use * to indicate a convolution defined as (/ * g)(x) — 



dp = -i[H FN {u s ),p]+~j / dxV[M v (x)]p + 7 / dxC[M v {x))p + ^7 / dxH[M v (x)]p <q{x) 
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J^ aQ dyf(y)g(x — y). This kernel function depends on 
the resolution length scale nu — x±/koXQ, where x± is 
the size of the condensate in the tight trapping direc- 
tions, xq = \JhJrruxi for loose trapping frequency w, and 
ko = 27r/A is the wavenumber of the probe light. Typ- 
ically, nu <C 1. The kernel function in Fourier space is 
given by fJ, v (k) = exp (— i^/c 4 ). 

In this system, the measurement signal provides a 
resolution-limited density measurement of the BEC. Un- 
like the cavity-mediated measurement, in no limit does 
the signal provide a position measurement. We investi- 
gate two values for v. v = 10 and v — 0.1, which provide 
'coarse' and 'fine' measurements of the BEC density, re- 
spectively. 

We examine the evolution of Eq. ^ after application 
of the Hartree-Fock method [12] : 

d t a H {x) = (-i(h F (x,u s ) + UN\a H (x)\ 2 ) 
+2-y(a u * /*„ * |a_H-| 2 )(;c) 
+Vl(Vv*V){vj)aLH(x), (7) 

where the Hartree-Fock order parameter an (x) is nor- 
malised to 1, so we must include a factor of N in the non- 
linearity term, accounting for the N atoms of the conden- 
sate. The NPW method converts Eq. (|Fj} to [131 |2"TI 151] : 

d t ajs/{x) = -i(h F (x,u s ) + U\a H (x)\ 2 

+Vl(fJ'v * T] f )(x)ja N (x), 

d t w = J dx (^({M^x)) 2 -2M v {x)A[M v {x)]) 
+2^fM„(x)n(x))w, (8) 

where M u (x) = (fi v * \a_\/\ 2 )(x). 

A comparison of the Hartree-Fock method to the NPW 
method for different strengths of v is shown in Fig. [3j 
(a). This analysis is performed without a quantum 
noise control. When v — 10 we see the Hartree-Fock 
method agrees with NPW, and both predict BEC cool- 
ing. This shows that for low resolution density mea- 
surements, higher-order modes are weakly excited such 
that feedback to the position mode provides sufficient 
damping, and the energy of the BEC stays low. Unfor- 
tunately, reaching this parameter limit is difficult with 
current traps [12]. When the BEC is probed at a finer 
spatial resolution (v = 0.1), the Hartree-Fock solution 
diverges from the NPW solution, which shows that the 
quantum noise of the spontaneous scattering events ex- 
cites higher-order modes of the BEC that are not cooled 
by the position control. This results in net heating. As 
with the cavity-mediated measurement, we aim to intro- 
duce controls for higher-order modes to compensate for 
the heating of the BEC. For this multi-channel measure- 
ment we shall require multiple channels of feedback, so 



we introduce an arbitrary time-varying potential for the 
atomic field. There are several proposals for implement- 
ing such potentials experimentally [35H39] . The single- 
particle Hamiltonian h F M{x,u s ,u^) describing this new 
feedback is: 

h F M{x,u s ,u^) = h F (x,u s ) + — 

x(n v *(j! u *Im[(fti,')})(x), (9) 

where is the strength of the quantum noise control for 
the direct measurement. 

The novel quantum noise control defined in Eq. ^ 
is simulated in Fig. (31(b). In part (b,i) we can see the 
novel control completely cancels the heating caused by 
the strong direct measurement of the BEC. The heating 
produced in a continuum of modes is efficiently damped 
by the field of controls. In parts (b,ii) and (b,iii) we 
compare the waveform of the BEC without and with the 
quantum noise control respectively. In (b,ii) we see the 
quantum noise causes the BEC to spread rapidly. In 
contrast, (b,iii) shows that the quantum noise control 
stops this distortion. The bulk of the BEC remains close 
to the centre of the trap. 

In both the direct measurement and the cavity- 
mediated measurement we have shown cooling of a con- 
densate is possible for a range of parameters with an ap- 
propriate choice of control. Experimentally, it is tempt- 
ing to try and cool a BEC in the 'weak' probing limit 
(£ = 0.1 for the cavity mediated and v = 10 for phase- 
contrast). In this case only the bulk feedback mecha- 
nism, governed by parameter u s and defined in Eq. ([I]), 
is required to achieve net cooling. However, operating 
in this parameter regime requires relatively small tightly 
trapped condensates [TT] [55] , although this may change 
with the development of trapping technologies. If the 
target condensate does not satisfy these requirements, 
more complex measurement and feedback is required. 
For moderately sized condensates a cavity-mediated mea- 
surement with an extra control appears promising. The 
shape of the quantum noise control suggested in Eq. ^ 
must match the measurement, and it might be possible 
to implement this feedback using the probe beam itself. 
For very large loosely trapped condensates, it is likely 
phase-contrast imaging will be the only option, in which 
case full spatial control of the trapping potential will be 
required. 

This letter has investigated the previously unexplored 
full quantum-field dynamics of a feedback-cooled multi- 
mode BEC undergoing (a) a cavity-mediated continuous 
cosine measurement, and (b) direct continuous density 
measurement via phase-contrast imaging. Full quantum- 
field simulations, performed with the recently developed 
NPW method, have revealed regimes with additional 
measurement-induced heating that is not included in 
single-mode and semiclassical models, and furthermore 
cannot be counteracted with simple linear and quadratic 
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FIG. 3: (a) Comparison of the Hartree-Fock method versus NPW for simulation of a cooled BEC under direct measurement 
for different v strengths. The NPW representation was integrated using Eq. Q, averaged over 10 paths, and plotted in red 
with circles. The Hartree-Fock representation was integrated using Eq. averaged over 100 paths, and plotted in purple 
with stars. A coarse resolution measurement, v = 10, is plotted in (a,i) while fine resolution measurement, v = 0.1, is plotted 
in (a,ii). Both plots have "f = 1, N — 100, U/N = 3 and u 3 — 1. The Hartree-Fock method only agrees with NPW when 
v = 10, it fails to predict the additional heating present in the system when ^ = 0.1 (b) A comparison of a controlled BEC 
under direct measurement without and with the quantum noise control averaged over 10 paths. (b,i) shows a comparison of 
the per particle energy of the system, while the density fluctuations in the system without and with control are shown in (b,ii) 
and (b,iii) respectively. In (b,i): integration of the NPW representation Eq. Q without quantum noise control, — 0, is 
plotted with red with circles; integration with the addition of the quantum noise control, u M = 5, is plotted in orange with 
with diamonds. (b,ii) and (b,iii) are contour plots of the BEC's density, all lines indicate equidensity parts of the BEC, the 
colour range of the lines range from (light) blue hues to (dark) red hues to indicate low density to high density respectively. 
All integrations were performed with u a — 1, u — 0.1, U/N = 3, 7 = 1 and N — 100. 



feedback controls. Fortunately, we have developed more 
sophisticated quantum-noise controls for both measure- 
ment schemes that cancel this additional heating, and 
shown that this allows for successful feedback-cooling of 
the condensate. This has been an important demon- 
stration of the necessity of full quantum-field simula- 
tions, and in particular the NPW method, to the suc- 
cessful design and implementation of measurement-based 
feedback-control schemes for multi-mode quantum sys- 
tems. 
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